Logarithmic Clustering in Submonolayer Epitaxial Growth 
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We investigate submonolayer epitaxial growth with a fixed monomer flux and irreversible aggre- 
gation of adatom islands due to their effective diffusion. When the diffusivity Dk of an island of 
mass k is proportional to k~'^, a Smoluchowski rate equation approach predicts steady behavior 
for < < 1, with the concentration Ck of islands of mass k varying as k~^^~^^^^. For /i > 1, 
continuous evolution occurs in which Ck{t) ~ {lnt)~^^''~^'^'^^^ , while the total island density increases 
as N{t) ~ (Int)'^^^. Monte Carlo simulations support these predictions. 
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Epitaxial thin film growth involves deposition of atoms 
onto a substrate and diffusion of these adatoms, leading 
to their aggregation into islands of ever-increasing size 
The resulting island morphology and mass distri- 
bution depends intimately on the diffusion processes of 
the adatoms. While this connection has long been rec- 
ognized [Q, a complete understanding of this evolution 
is still lacking. For a variety of atomic transport mech- 
anisms, there is a power-law dependence of the effective 
island diffusivity Dk on its mass k, Dk oc fc"'' [§^7 with 
fi typically in the range (1/2,3/2). Starting with this 
observation, various approaches have suggested that the 
growth of islands due to their clustering is a power law 
in time 

In this Letter, we provide a comprehensive account 
for the evolution of the island size distribution in the 
submonolayer regime by solving the Smoluchowski rate 
equations. For mobility exponent < fi < I, a steady 
state arises in which the concentration of islands of mass 
k is given by Ck oc with t — {3 — A*)/2. For 

all /i > 1, logarithmic island evolution occurs in which 
their total density grows as (Int)^/^ while Ck{t) varies as 
(Int)"'^^'^"^-''^/^. More generally, our approach applies to 
any epitaxial system in which the diffusivity of an island 
vanishes more rapidly than inversely with its mass. 

In the diffusion-controlled limit, the aggregation rate 
Kij of an i-mer and a j-mev is given by the Smoluchowski 



formula Kij ~ {Di + Dj){Ri 



ra. Here R, is 



the linear size of an i-mer (island of mass i), which is as- 
sumed to be compact, and d is the spatial dimensionality 
of the substrate. This Smoluchowski formula is applica- 
ble in d > 2, while in the physically relevant case of two 
dimensions the reaction rate depends logarithmically on 
the island size [|l^. For both simplicity and because lit- 
tle quantitative information is lost, we shall ignore these 
logarithmic factors. This is equivalent to treating the 
islands as point-like throughout their evolution. 

With Dk oc k~^^, an appropriate choice of time units, 
and the neglect of logarithmic corrections, the reaction 
rate in two dimensions becomes 



(1) 



and the corresponding Smoluchowski rate equations for 
the evolution of the concentrations of fc-mers in the pres- 
ence of a steady monomer flux are 
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These rate equations represent a mean-field approxima- 
tion in which spatial fluctuations are neglected, and also 
a low-coverage approximation, since only binary interac- 
tions are treated. 

Let us first consider the behavior in the steady state 
regime. To solve the rate equations in this case, we in- 
troduce two generating functions 



oc 



C{z) = } ^ CkZ- 



00 
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(3) 



Multiplying Eq. (0) by z^, and summing over all k, gives 



C^iz)Ciz) - C^iz)N - Ciz)N^ + Fz = 0. 



(4) 



Here = C(z = 1) = ^ cj, is the total island density 
andiV^ = C^(2-l)-^fc-'^Cfc. 

We now assume a power law asymptotic behavior for 
the steady state concentration, 



C 



(5) 



as fc ^ cxD. For this power law to hold for all fc, we re- 
quire T > 1, so that J2 converges; this leads to the 
condition fi < 1 for the mobility exponent as shown be- 
low. From basic Tauberian theorems , the asymptotic 
form for Cfc in Eq. (|^) induces the following power-law 
singularities in the generating functions as 2; ^ 1 



C{z) =N + CTil - r)(l - zy-^ + 



Cf,{z) ^ + cr{i 



fi)il-z) 



(6) 



The leading constant factor in each line is finite and coin- 
cides with the definition given in Eq. (|) if the exponent 



1 



of the second term is positive. Otherwise, the constant 
factor vanishes and the generating function has a power- 
law divergence as z — > 1. Substituting these expansions 
into Eq. (jj) and matching the leading behavior in (1 — z) 
leads to the decay exponent t = (3 — /i)/2. The condition 
for a steady state to occur, t > 1, thus imposes an upper 
bound on the mobility exponent, /i < 1. From matching 
the leading behavior in (1 — z), the constant C may also 
be determined, from which the island mass distribution 
in the steady-state regime < /i < 1 is 



Cfc 



47r 



(7) 



It is important to note that this mass distribution holds 
only up to a mass cutoff kc{t) ^ whose value is deter- 
mined by requiring that the total mass in the system due 
to the steady input is proportional to t - islands of mass 
greater than kc{t) have not yet formed. Therefore 

M(t) = f2 kck{t) ^ fc^'^"'^/" ~ ^ (8) 



fc=i 



fc=i 



which gives the mass cutoff exponent 
C(Ai) =2/(At + l). 



(9) 



We now investigate the asymptotic behavior of the is- 
land mass distribution for fi > 1. In the extreme case 
of n = oo, i.e., diffusing monomers and immobile is- 
lands, the island density grows as a power law in time. 
Nit) ~ (3^2^)1/3 We shall argue that contin- 

uously evolving behavior occurs for all /i > 1, but with 
anomalously slow logarithmic kinetics. When /i is strictly 
greater than unity but still finite, we find 



N{t) 



sin(7r//x) 



InT 



m/2 



(10) 



with T = t\/F, while the concentration of /c-mers decays 
in time as 



(11) 



It is remarkable that such logarithmic dependences, a 
feature which generally signals marginal behavior, oc- 
curs in the entire regime 1 < ^ < oo. In the borderline 
case of /i = 1, even more unusual behavior arises with 
N{t) - y/lnT/ In(lnr). 

Our argument leading to Eqs. and ([ll| ) is based on 
a quasi-static approximation, in which the time deriva- 
tive in Eq. (^) is neglected. Indeed, the logarithmic be- 
havior in Eqs. ( p^ and ( pi] ) immediately implies that the 
temporal derivatives in the Smoluchowski rate equations 
are asymptotically negligible. Within this quasi-static 
framework, Eqs. (0) become 



ci [N + N^) 



E 

i+j—k 



{i-i^ + r^") c^c, - Ck (k-m + N^) . (12) 



Further, by summing Eqs. (^2|) over all k, the total island 
density in the quasi-static limit obeys 



= 1- NN„, 



(13) 



In Eqs. (^2|) and ([l^) we have set = 1 by a rescaling of 
units. Eq. (|l^) immediately gives = N^'^, and then 
from the first of Eqs. (p2|), ci ~ 1/7V. The remainder 
of Eqs. ( p2|) may then be solved recursively. By writ- 
ing the first few of these equations, it is evident that the 
dominant contribution to is the term in the quadratic 
product which is proportional to ciCk-i- If we keep only 
this contribution, the resulting recursion may be solved 
straightforwardly to yield 



fc-i 



Ck 



k k-1 

i=2 j=i 



+ f 



(14) 



Since the factors Bj < 1 for < TV^, while Bj 1 
for j^^ ^ TV^, this implies that Ck is a rapidly decreasing 
function of k for k ^ A^^ and then becomes constant for 
larger k. 

To compute Ck, first note that for /i > 1 the product 
Ylj bj converges, so that it may treated as constant. We 
then write the second product as the exponential of a 
sum and take the continuum limit. This leads to 



Cfc ~ exp 



TV 



exp 



J=2 



-iV^/'' / ln(l + iy-^)dw 
/o 



(15) 



where w = j/N'^/'^ and x = k/N'^^^. This form has 
two slightly different asymptotic behaviors depending on 
whether ^ is strictly greater than or equals 1. For /i > 1, 
the monotonically increasing integral in Eq. (Hq) con- 
verges as X ^ oo. Thus Ck decreases as a function of 
k until a threshold value fcth — N'^^'^, beyond which Ck 
remains constant with a value determined by taking the 
upper limit of the integral as infinite. Hence 



Cth 



(16) 



with = Jq°° ln(l-|-w;~'^) dw = tt/ sin(7r//x). Physically, 
we expect this constancy in Ck to persist until k reaches 
the cutoff kc ^ t'' . 

To check this result, we performed numerical simu- 
lations in the mean field limit of submonolayer epitax- 
ial growth. In the simulation, an island of mass fc. 



2 



which remains point-Hke throughout the aggregation pro- 
cess, moves equiprobably to any site with a probabil- 
ity proportional to fc"'', as mandated by the power-law 
mass-dependent island difFusivity. There is also a steady 
monomer flux entering the system. As shown in Fig. 1, 
Ck{t) is nearly constant in k over a substantial range as 
predicted by our theory. However, when k kc, there is 
a peak in Ck which is not accounted for in our quasi-static 
description. 




FIG. 1. Ck{t) versus on a double logarithmic scale at 
t 22000 for fi = 1.5. The data is based on 5000 realiza- 
tions of an initially empty system with F = 0.05. 
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FIG. 2. N(t) versus ln{tf^'^) for = 1.2 (o), 1.4 (o), 1.5 
(A), and 2.0 (v)- The data is based on 1000 realizations. 
Also shown are power law fits to the data in the range t > 3. 
This gives, for the exponent of Int, 0.38, 0.55, 0.64 and 0.76, 
respectively, for ^ = 1.2, 1.4, 1.5, and 2.0. 

The time dependence of the total island concentra- 
tion may now be determined by using the fact that the 



mass density in the system grows linearly with time, 
M{t) — J2k>i ^^'^k{t) = t. The dominant contribution 
to this sum is from the plateau region fcth < k < kc, 
where is approximately constant in fc, Ck ~ Cth- Thus 



M{t) 



Cth ^ k . 

k=kth 



(17) 



Here we also use the fact that kth grows only logarith- 
mically in time (see below), so that the lower limit can 
be taken to be zero. Furthermore, kc{t) ~ grows 
more rapidly than t^^^. Indeed, since the cutoff exponent 
C(w) = 2/{ii + l) for At < 1 (Eq. (|)) and ({fi = oo) = 2/3 
|12|-p^, we anticipate that 2/3 < C(//) < 1 when 1 < /i < 
oo. In fact, C = 1 in this range of mobility exponent, as 
we will show below. Using C = 1, Eq. ([T^), together with 
Cth exp(— A)jA^^/^) and kc ^ i^, we immediately ob- 
tain Eq. (pl)|). Our data for N{t) as a function of time 
(Fig. 2) is qualitatively consistent with N{t) growing as 
a power of Ini, but with a smaller exponent than fi/2. 
Note finally that fcth ^ N'^/'^ which is proportional to 
Int. 

To determine Ck{t) for fc^ ^ N"^, we now use the ap- 
proximation Bj sa jf^/N^ in Eq. (p^ to give 



Cfc 



(18) 



which directly leads to Eq. ([Til). Finally, the time de- 
pendence of Cth (t) and kc (t) may be determined from the 
sum rules (0) and J2 Ck ~ cthfcc ~ N. These two rela- 
tions give Cth ~ N'^ /t and kc ^ t/N. Thus in the plateau 
region fcth < k < kc, 



Ck(t) 



{inty 



kc{t) 



(19) 



It is important to note that our approach applies to 
any mass-dependent island diffusivity which decays faster 
than its inverse mass. For this general situation, the ana- 
log of Eq. (0) is Ck ^ + Dj)il + N^Dj)-\ 

For example, for Dk ~ e"^"*^, a case that wa s con sidered 
numerically in |l^, we obtain N{t) ~ exp(-\/a hit). This 
unusual growth - faster than any power of logarithm but 
slower than any power law - would be difficulty to ob- 
serve numerically. 

In the specific case /i = 1, subtler nested logarith- 
mic behavior arises, as reflected by the additional sin- 
gularity in Eq. ( p^ as /i 1. First, the product 

11^=1 = rii^i^ll + J"^) in Eq. (^ now equals k. Sec- 
ond, the term C2Cfc_2 also contributes to the asymptotic 
behavior. Third, and most importantly, the integral in 
Eq. (|l^) diverges at the upper limit. Due to the second 
attribute, the recursion relation for becomes 



Cfc 1_ 

k 



k/N^ 



k/N^ 



Cfc-l 

fc- 1 



Cfe-2 



1 

2 iV2- 



We seek a solution for Cfc in the form of Eq. 
we write 



(20) 
M). Thus 



3 



Cfe 



(21) 



i=2 



where the factor Ck accounts for the additional term in 
Eq. (|^). Substituting into Eq. (|2^) gives 



Ck — Cfc-i + C, 



fe-2 



k-1 



1 



(22) 



These coefficients are slowly varying in k when /c 3> 1 and 
we may treat k as continuous in this asymptotic regime. 
Eq. (p2) then becomes a differential equation whose so- 
lution is Cfc ~ fce^ (with x = k/N'^). Consequently, 



Cfe ~ exp 



ln(l 



dw 



(23) 



Thus for /i = 1, Cfe decreases rapidly in A: for A: ^ iV^ 
and then increases until k reaches k^. The island mass 
distribution attains a minimum at xth = N"^ whose value 
is 



cth exp 



IniV^l 



(24) 



Paralleling the analysis of the case /i > 1, the total 
island concentration is 



N{t) 



' hit 
In(lni)' 



(25) 



while Eq. ( pT| ) together with Ck ^ k implies that the 
concentration of islands for mass k ^ iV^ is 



Cfe 



(fc + 1)! 



(fc + 1)! 



In(lni) 
hit 



fe-l/2 



(26) 



Finally from Eq. (|l9|), the time dependence of Cfc(i) and 
fcc(i) is given by 



Cfc(i) 



1 \nt 



t In(lnt)' 



fcc(t) t 



In(lni) 
Int 



(27) 



Our results should be generally applicable to real epi- 
taxial systems in the submonolayer regime. This regime 
requires -Ft <C 1, while the asymptotic predictions of our 
Smoluchowski theory apply for t^/F ^ 1. Consequently, 
our results should be valid for t F^^. Since 

the dimensionless flux F is small in typical epitaxy exper- 
iments, the time range over which our theory will apply is 
correspondingly large. Notice that the maximum island 
density attained at the end of the submonolayer regime 
imax - F-^ scales with flux as iV^ax - F^/'^[hi{l/ F)Y/'^ . 
In fact, for all systems with the diffusivity of an island 
decaying more rapidly than its inverse mass, our ap- 
proach leads to iVmax universally being proportional to 
F^/"^ times a subdominant model-dependent factor. 

In conclusion, we determined the kinetics of island- 
ing in submonolayer epitaxial growth, in which adatom 



hopping induces a power-law mass-dependent island dif- 
fusion, with Dk oc k~^^. This leads to the reaction 
rate Kij oc {^^^^ + j^^) between two islands of mass 
i and j in the Smoluchowski rate equation. A steady 
state arises for mobility exponent < M < 1, in which 
the island concentration varies as Ck ^ fc^*^"^^^^/^ for 
k < kc oc with C{^i) = 2/(1 + A*)- Strikingly, log- 
arithmic time dependence arises for all 1 < /Lt < oo, 
a feature suggestive of marginal behavior over this en- 
tire range. In this regime, the total island density N{t) 
grows as (Int)'^/^, while the density of islands of mass fc 
is Ck{t) cx (lnt)-(2fc-i)A'/2^ for < fcth oc \nt, Ck{t) inde- 
pendent of fc, with Ck{t) ^ t~^{h\t)^^, for fcth < k < kc, 
and Ck vanishingly small for k > kc- For /i = 1, even 
more unusual nested logarithmic behavior occurs. 
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